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摘 要 :为 了 进一步 提高 求解 Volterra 型 积分 微分 的 数值 精度 ,针对 一 种 变 系数 Volterra 型 积分 微分 方程 ,提出 了 2 种 

Legendre 谱 Galerkin 数值 积分 法 。 采 用 Galerkin Legendre 数值 积分 对 Volterra 型 积分 微分 方程 的 积分 项 进行 预 处 理 , 对 
< 其 构造 Legendre tau 格式 ,同时 用 Chebyshev-Gauss-Lobatto 配置 点 对 变 系数 和 积分 项 部 分 进行 计算 ,并 通过 对 方程 的 定义 
bent. 提出 了 一 种 多 区 间 Legendre 谱 Galerkin 数值 积分 法 。 该 方法 的 格式 对 于 奇数 阶 模型 具有 对 称 结构 。 此 
人 外 ,通过 引入 Volterra 型 积分 微分 方程 的 最 小 二 乘 函 数 ,构造 了 Legendre iif Galerkin 最 小 二 乘 数值 积分 法 。 该 方法 对 应 的 
e 数 方程 系数 矩阵 是 对 称 正 定 的 。 数 值 算 例 验 证 了 这 2 种 Legendre 谱 Galerkin 数值 积分 方法 的 高 阶 精度 和 有 效 性 。 
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"Abstract. In order to further improve the numerical accuracy of solving Volterra integro-differential. two kinds of Legendre 
spectral Galerkin numerical integration methods are investigated for the Volterra-type integro-differential equation with vari- 
able coefficients. Firstly, the Galerkin Legendre numerical integration is applied to deal with the integral term of the Volter- 
ra-type integro-differential equations. Secondly. the Legendre tau scheme is developed for the Volterra-type integral-differ- 
ential equations with variable coefficient, and the Chebyshev-Gauss-Lobatto collocation point is used to the calculation of the 


variable coefficient and integral term. Finally. by decomposing the definition interval of the function. the multi-interval Leg- 


endre spectral Galerkin numerical integration method is also designed. Its scheme of the proposed method has symmetric 
structure for odd-order model. In addition, by introducing the least squares function of the Volterra type integro-differential 


equation. the Legendre spectral Galerkin least-squares numerical integration method of is constructed. The corresponding 


coefficient matrix of the algebraic equation is symmetric positive. Some numerical examples are given to test the high-order 
accuracy and the effectiveness of our methods. 
Key words; integral differential equations; numerical integration; Chebyshev-Gauss-Lobatto interpolation; least square meth- 


od; Legendre Galerkin 


收 稿 日 期 : 2022-03-31 

基金 项 目 : 国家 自然 科学 基金 (12161025) ;广西 自动 检测 技术 与 仪器 重点 实验 室 基金 (YQ22106) ;广西 科技 基地 和 人 才 专 项 ( 桂 科 AD18281025) ; 
桂林 电子 科技 大 学 研究 生 教育 创新 计划 (2020YCXS086) 

通信 作者 : 覃 永 辉 (1985 一 ), 男 ,副教授 ,博士 ,研究 方向 为 偏 微分 方程 数值 解法 。E-mail:yonghuil1676@163. com 

引文 格式 : 范 友 康 , KLAR. KE. 变 系 数 Volterra 型 积分 微分 方程 的 2 种 Legendre 谱 Galerkin 数值 积分 方法 [J]. 桂林 电子 科技 大 学 学 报 ， 
2023 ,43(2) :158-164. 


第 2 期 


Volterra 型 积分 微分 方程 作为 带 记忆 性 质 的 数 
学 模型 ,通常 出 现在 超 流 体 问题 "、 金 融 数 学 问题 
4 Turbulent 扩散 问题 汪 等 领域 。 通 常 来 说 ,由 于 
Volterra 型 积分 微分 方程 具有 记忆 性 质 , 它 的 精确 解 
很 难得 到 ,上 且 也 为 构造 其 数值 解法 带 来 许多 麻烦 。 
为 了 推广 Volterra 型 积分 微分 方程 在 实际 问题 的 应 
用 ,相关 学 者 进行 了 研究 5 。 另 外 ,对 于 变 系数 Volt- 
erra 积分 微分 方程 , 当 采 用 谱 配 置 方 法 允 近 时 ,得 到 的 
和 矩阵 往往 是 稠密 的 ,因此 增加 了 求解 的 难度 中 。 

考虑 如 下 积分 微分 方程 

put =a (x)ulx) + f(x) + 


4 | KGa ds, seg] 0 


u (0) 一 Wo。 

E% D :—(Gslj0osisszmz T) 假设 f(x)、 
D 为 已 知 的 光滑 函数 , K (x ,s) 为 核 函 数 。 对 于 
ARA Voltrra 积分 微分 方程 (1) ,近年 来 ,各 种 数值 
4 也 得 到 广泛 关注 与 研究 , A BO EA, ARR 

ge-Kutta 1" * hp- 连 续 Petrov-Galerkin 有 限 元 
ILE 75 157 Legendre Galerkin 77r 35^" , if 
rkin 多 项 式 法 0 、 谱 Petrov-Galerkin 7; z^, 
ing Íg/]N 3e; Ag. 
许多 学 者 对 谱 方 法 求解 积分 方程 进行 了 研究 。 
tfi lang" 人 研究 了 Volterra 型 积分 微分 方程 的 谱 
Legendre 配置 法 ,从 理论 上 给 出 了 谱 精 度 收敛 的 证 
Wee Li 等 中 给 出 了 线性 重心 有 理 法 求解 Volterra 积 
4587; E. Gu 介绍 了 一 种 多 步 Chebyshev TÉ 
配置 法 。Wang 等 中 研究 了 Volterra 积分 微分 方程 
的 Ship-Legendre Jacobi 谱 配置 法 ,得 到 了 关于 光滑 核 
和 能 奇异 核 的 一 些 近 似 结果 。 

Tang 等 "研究 了 具有 周期 边界 条 件 的 抛物 型 
方程 的 单 区 间 和 多 区 间 Legendre tau 谱 方 法 。 
Zhang 等 ”针对 一 阶 常 微分 方程 的 初 值 问题 提出 了 
Legendre tau 方法 。Qin 等 5 研究 了 发 展 方程 的 
Legendre tau Galerkin 谱 配 置 法 。Karimi 等 5 提出 
了 一 种 最 小 二 乘法 求解 线性 Fredholm 积分 方程 。 
胡 玉 巧 等 外 针对 一 种 抛物 型 方程 研究 了 Legendre 
Galerkin 谱 配 置 最 小 二 乘法 。Maleknejad 457" Ji tH 
了 一 种 块 脉冲 函数 运算 矩阵 法 , 它 利用 块 脉冲 函数 及 
其 积分 的 随机 运算 矩阵 ,可 将 Ito-Volterra 积分 方程 
简化 为 线性 下 三 角 方 程 组 ,但 得 到 的 结果 精度 不 是 很 
高 。 在 文献 L[25] 的 基础 上 ,Ahmadinia 等 .站 提出 了 
一 种 基于 最 小 二 乘法 和 块 脉冲 函数 的 求解 Ito- Volt- 
erra 积分 方程 的 计算 方法 ,其 精度 优 于 块 脉冲 函数 运 
算 抢 阵 法 。 本 文 的 方法 对 积分 微分 方程 中 的 积分 项 


范 友 康 等 : 变 系数 Volterra 型 积分 微分 方程 的 2 种 Legendre 15 Galerkin 数 值 积分 方法 praes 159 


CIHITIGATV c T FRATI 
部 分 先 采 用 Galerkin Chebyshev-Gauss-Lobatto 或 
Legendre-Gauss-Lobatto 数值 积分 方法 进行 离散 , 然 
后 给 出 Legendre tau 格式 ,采用 Chebyshev 对 变 系 
数 部 分 和 源 项 进行 插值 处 理 。 
1 符号 和 格式 
对 任意 的 > > 0, 设 H'A) 为 经 典 Sobolev 空 
间 。 令 PT) 为 次 数 不 超 过 N 的 多 项 式 集合 。I% 
为 Chebyshev-Gauss-Lobatto 插值 算 子 IS :C (1) 一 
Pa), 使 得 
TAu (xz ) 一 xxz)，0 委 ;过 N。 
S [&1C GOD — PNGT) 为 Legendre-Gauss-Lobatto 
插值 算 子 , 使 得 
I&uCx) =u), 0Zi<N, 
其 中 , {xf} Get) 4:323 Chebyshev-Gauss-Labotto 
点 、Legendre-Gauss-Labotto 点 。 


1.1 Legendre i Galerkin 数值 积分 法 


考虑 式 (1) 的 弱 格 式 : 求 x € H'A), 使 得 
MTM 
T Ksuls ds) + (fs0), (2) 
u (0) = Uo» Vv c L Oi 
因此 , 式 (2) 的 Legendre 谱 Galerkin 数值 积分 法 为 : 
对 于 任意 v € Pya Cl), 求 UN c Pn(1) ,使 得 
me =(IŞ[aun]; v) + Afv) + 


1 aS IG Dun Gds s0), (3) 


u(0) 一 xu 。 
考虑 多 区 间 Legendre 谱 Galerkin 数值 积分 格 
式 。 首 先 , 将 区 间 工 =L0,T] 划分 为 一 些 子 区 间 
Q— xm, eM €um T, 
即 1-0 Gam]. P—1.2.5M, 

u (x) =u(x) sete ugs d 过 iM。 

邻 T,= 二 (xz; szi] i = 二 1,2,…,M。 设 Pn, (1;) 
为 区 间 I; 上 次 数 不 超 过 NN; 的 多 项 式 集合 。 定 义 I, 
上 对 应 的 Chebyshev 插值 算 子 IL COD 一 
Py dD X 

AEn u) |i, Ep Su lr, G£,),j 90, s N; 
zi € Tale 0e M, 

Legendre 插值 算 子 It, COO Py OON 
Ainu) | 1 (zk;) 5u | 1 (215743 =0,1,2,.,N;, 
zj EL,l1<i<M., 
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在 每 个 区 间 采 用 Legendre 谱 Galerkin 数值 积 
分 求解 , 即 式 (1) 的 多 区 间 Legendre Galerkin 数值 积 


分 格式 为 : 

对 于 任意 v € Pua da), SR uh € Py, da) ,使 
4H 
F 


uk, o) — OE, Laut, o) + AEn, fov) t 
k—1 
C 
> ats], 


AEn | Tin, CCG i, GI ds D. 


FV K(x puk, (dy. v) + 


un =u (0) un SuN Gia). 


(4) 
1. 2 算法 实施 
TORREK UU 
C GO =1, $i G) =Lro (£) Lq14G);. 
gi) GO =Lra G)— Less. 2k <N, 


jp. L4 GO 为 位 移 Legendre 多 项 式 ; 1 为 多 项 式 


N, 
= uN, G) = Opi zr) Huk, o 


ee i 一 1 
JORIE RE S 和 质量 矩阵 M : 
Y S —G;) —G'jP Ly, 
~ M —OM;)-—(C3Sj Lr). 
= } ata 
Df =r | a Ee. =s (x idu. 
区 | 2 2 
(99) 
CERU) P, Wv =Lr,: (z), L= i; Zata N—1; 
lll (4) 的 离散 系统 为 
($5— M* — A“*)u’ =F, u = (kih, th) 9 
其 中 ， 


Is [^ Thu, QCG Du, GDds]— 


ls 5» Y KGosG 09; GG sv; Doy s 
i=l j=0 


ON N 
tm E» S KG str D$, (sG sr; Po, iu) , 


A ie 
F= IŞ fL) p» übs | s Kou ds Ln), 
oren- | 

乘 数 值 积 分 法 


1.3 Legendre 谱 Galerkin 最 小 二 


对 于 任意 的 wE PG 式 (1) 的 最 小 二 
可 定义 为 


s G)4 


Fw; Ð= |F 本 
| KG DoGods l2, (5) 
相应 的 最 小 值 问 题 为 : 求 vE PaO), 使 得 
FCu; f) p ci 
因此 , 式 (5) 的 变 分 形式 为 : 
Ru € PO), 使 得 
(2,u,9,v) — (aCx)uCx).9,v) — 
È Keu Oda (3,0. Vv € PC. 
(6) 
一 步 ,给 出 式 (6) 的 Legendre Galerkin 最 小 二 
乘 数值 积分 格式 : 求 ww € PaO), 使 得 
(Qux 2,0) — (OISlaGOuxs (GO 9,0) — 
as]. I&K Gros)ugGOdss9,9) — US f 8,05 
Vv € Pll). C) 
类 似 地 ,基于 上 述 的 区 间 划 分 , 式 (7) 的 多 区 间 
Legendre Galerkin 最 小 二 乘 数值 积分 格式 为 : 
求 uN, € Py GO, 使 得 
(QuN, ,9.0) = tu, Laui, 00) H Ain, ao) 十 


£1 
P? bs]. Iis, ts Du, (QD dg.2,v) 十 
i=1 i 


AEn [ In, K @s)uk, )ds,30), 
T, bi T 


un =u (0) uh =uh (z); Yv E Py, CHP 


(8) 
格式 (7) 的 实施 过 程 : 令 基 函数 为 如 下 形式 "1 

eh 

pi(z) =E a HL E), 

P(x) —L,G) — LG») sk zm 2 
引进 和 矩阵 

W —OY;)-—(2,9;.299;); 

^ 也 — 9, ;过 二 0,1,.…,N， 并 代入 式 (7), 可 得 


Qy —IDu =F, 
其 中 : 
II — (II) — (I Sa G9, ,991) 十 


aS KG soe Gods ,9.90); 
F —USf.2,9) 
2 数值 算 例 
分 别 给 出 一 些 数值 算 例 验 证 所 构造 的 2 种 方法 
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的 高 阶 精 度 与 有 效 性 。 利 用 2 种 数值 格式 解 出 它们 
的 数值 解 ,并 与 精确 解 进行 比较 , 求 出 工 ”- 范 数 误 
差 。 将 2 种 方法 的 数值 实验 得 到 的 数值 误差 与 文献 
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Ti HT! 


误差 的 图 像 。 从 表 1 和 图 2 可 看 出 , 当 m=3 时 , 式 
4) 的 精度 略 优 于 MCSC 方法 的 精度 ;mx —1 时 , 式 
4) 的 误差 与 MCSC 方法 类 似 。 因 此 ,可 以 断定 式 


PR ^ 


L18j] 中 的 误差 进行 比较 。 定 义 (4) 对 于 Volterra 积分 微分 方程 是 有 高 阶 精 度 的 。 
h = max C; — x; 4)2/2:i —0,1,*. M), 表 1 N,=N=3 时 , 例 2 的 L™ - 范 数 误 差 
例 1 考虑 Volterra 积分 微分 方程 
9020 C€H5) 方法 m i 
aCx)-coszc; KZ 一 了 Te 5 . (9) 50 100 200 
de i o 1 4.22X10 ^  1.49x10 ' 5.27X10 ^? 
应 的 精确 解 为 u(x) = NE 3 — 2.01X107  1.83X107  1.65X10^? 
本 例 验证 Legendre Galerkin 数值 积分 法 的 有 效 1.59X10™*  5.63X10 ^ 1.99X10™ 
j n > DE , Hr 的 M 格式 (4) 
E. 分 别 用 式 (3)、(4) 求 解 本 例 , 将 单 区 间 的 Leg 3.00 10^? 2. 65x 10^? 2.35x]0-? 
endre Galerkin 数值 积分 式 (3) 和 多 区 间 的 式 (4) 进 
行 比 较 , 区 间 划 分 个 数 M 取 1、2, 其 误差 如 图 1 所 - -O -- MCSC,n-1 
示 。 从 图 1 可 看 出 , 当 多 项 式 次 数 N 取 10 时 ,精度 NE rcoem 
T 107" ,这 表明 该 方法 具有 高 阶 谱 精 度 ; 当 取 相 LC E ZEE 
Mts IRU N M, M2 时 得 到 的 精度 优 于 M To M HEEEEE-EEEERET: 
T aps rats 
à ~N 
Bau. uus 
- - 米 - M=1 = XK x 
--Q- M2 USB e n 
` 50 100 150 200 
cA 10x 1/200 
AN 图 2 例 2 的 L”- 范 数 图 像 
Se 米 、 
一 例 3 考虑 具有 如 下 核 函数 和 精确 解 的 积分 微 
6 8 10 12 14 , 
N KGr,s)—-— pui "legit D, ax)-—1. 
图 1 例 1 的 工 ”- 范 数 图 像 zT (11) 


例 2 考虑 Volterra 积分 微分 方程 


Qu- 2m tl o QU sinx l gamn | 
2 ` 2(m +1) l 
| sin xu +Ë zu +Ë, ma uCs2ds. 
(10) 
其 精确 解 为 wx(z) 一 z ^. € (0,2]。 


对 于 例 2 ,文献 [18] 用 MCSC 方法 求解 Volterra 
积分 微分 方程 (10), 且 在 表 1 中 给 出 了 其 工 ”- 范 数 
误差 。 利 用 式 (4) 求 解 例 2 Pm =1,3 时 的 Volterra 
积分 微分 方程 (10) , 表 1 为 工 ”- 范 数 误差 。 为 了 方 
便 2 种 方法 进行 比较 ,依次 取 1/4& 二 50,100,200 和 多 
项 式 的 次 数 为 N; =N ==3。 图 2 为 它们 的 L”- 范 数 


初 值 条 件 为 u(0) —0, 

在 例 3 中 ,利用 多 区 间 Legendre 谱 Galerkin 数 
值 积分 法 CAO 处理 计 算 区 间 较 大 (T 从 50 增加 到 
103) 的 问题 , 取 N= 二 10, k=1/2. 对 应 的 LL”- 范 数 
误差 如 图 3 所 示 。 从 图 3 可 知 , 随 着 了 的 增 大 ,误差 
也 随 着 增 大 , 当 区 间 长 度 增 大 到 全 二 350 时 ,误差 稳 


定 在 10 ,因此 当选 择 合适 的 多 项 式 次 数 和 区 间 划 
分 ,该 方法 对 计算 区 间 较 大 的 问题 也 有 效 , 且 具有 高 


阶 谱 精 度 。 
例 4 考虑 如 下 积分 微分 方程 
Le sinx —ax)4 


(sin z2uCx) 十 
4 | e7«coas, x € (0.2], 


Lu =l 
(12) 
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假设 其 精确 解 为 u(x) 二。 


0 0 
| k =1/2, N =10 


"| 


y 


1 1 上 上 L J 
50 200 350 500 650 800 1000 


llu-u loo 


T 
3 例 3 的 L”- 范 数 误差 


2o. 45051 FH A JE. Chebyshev 谱 配 置 (MCSC) 方 
ks Art s 5] ,其 结果 如 文献 [18] 中 图 1 所 示 。 图 4 为 
As L”- 范 数 误 差 曲 线 。1/ 二 1 时 ,MCSC 对 应 
单 步 法 , 取 多 项 式 次 数 N —2. 4. 6. 8. 10 和 
EXT 4 可 知 ,2 种 方法 的 元” - 范 数 误差 相似 。 


- -© - MCSC 
- 3k - 格式 (7) 


4 例 4 关 于 式 (7) 和 MCSC 法 的 L”- 范 数 误 差 图 像 

例 4 还 考虑 了 MCSC 方法 取 不 同 区 间 划 分 1/= 
10,50,100,150 和 200 时 的 情况 , 且 多 项 式 的 次 数 为 
N: SN =3, 对 应 的 L”- 范 数 误 差 如 图 5 所 示 。 为 了 
表明 式 (8) 的 谱 精 度 ,在 计算 时 同样 取 1/4 — 10.50, 
100,150 和 200, 图 5 为 它 的 误差 。 该 方法 给 出 的 数 
值 结 果 略 优 于 文献 [18]。 

例 $ 考虑 Volterra 积分 微分 方程 


10 [ 
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分 别 利用 文献 [18j 中 的 MCSC 方法 和 多 区 间 
Legendre 谱 Galerkin 最 小 二 乘 数值 积分 法 (8) 求 解 
式 (13) 的 L”- 范 数 误 差 ,其 图 像 如 图 6 所 示 。 文 献 
[18] 中 ,配置 点 固定 数 为 11( 即 多 项 式 次 数 N =10), 
子 区 间 大 小 固定 为 = 二 1/2。 为 了 方便 比较 , 令 式 
(8) 计 算 时 多 项 式 次 数 固定 为 N 王 10( 配 置 点 数量 为 
1D. 3 —1/2, & T=2,5,10,15,20. MBI 6 可 看 
出 , 式 (8) 在 工大 于 10 时 的 精度 略 优 于 MCSC 方法 
的 精度 ,表明 该 方法 在 计算 区 间 略 大 的 问题 时 是 有 优 
势 的 。 


3 RE 


对 变 系数 Volterra 型 积分 微分 方程 ,分 别提 出 
了 Legendre 谱 Galerkin 数值 积分 方法 和 Legendre 
Galerkin 最 小 二 乘 方法 。 第 一 种 方法 采用 Legendre 
tau 方法 进行 离散 ,因此 ,对 奇数 阶 方程 其 格式 具有 
对 称 的 特点 ,从 而 方便 今后 格式 收敛 性 分 析 的 研究 。 
在 计算 中 ,采用 了 Chebyshev-Gauss-Lobatto 配置 点 
插值 ,使 得 格式 可 采用 FFT 算法 得 到 其 系数 ,所 以 提 
高 了 格式 计算 效率 。 对 区 间 较 大 的 模型 ,通过 合理 选 
择 区 间 剂 分 数量 和 增加 多 项 式 次 数 , 使 得 算法 在 实施 
中 具有 自 适 应 技术 的 优点 。 第 二 种 方法 ,基于 最 小 二 
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6 例 5 的 L”- 范 数 误 差 


乘 原理 ,人 研究 了 Legendre 谱 Galerkin 最 小 二 乘 数值 
积分 法 。 该 方法 对 应 的 离散 方程 具有 对 称 结构 ,因此 
便 邓 采用 迭代 方法 求解 。 数 值 算 例 验证 了 这 两 类 
Liiendre ilf Galerkin 方法 的 高 阶 精度 与 有 效 性 。 数 
值 冤 例 说 明了 这 两 类 方法 对 于 变 系数 Volterra 型 积 
f 身 分 方程 有 高 阶 精 度 。 此 外 ,将 该 方法 的 数值 实验 
得 到 的 数值 误差 与 文献 L18j 中 的 误差 进行 比较 , 结 
表明 该 方法 具有 一 定 的 优越 性 。 下 一 步 将 研究 将 该 
大 并 推 广 到 非 线 性 的 Volterra 的 积分 微分 方程 和 带 
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